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Oh' theory. Making use of Starobinsky's stochastic inflation we prove that deep infrared scalar 



We consider phase transitions on (eternal) de Sitter in an 0{N) symmetric scalar field 



modes cannot form a condensate - and hence they see an effective potential that allows 
no phase transition. We show that by proving convexity of the effective potential that 
governs deep infrared field fluctuations both at the origin as well as at arbitrary values 
of the field. Next, we present numerical plots of the scalar field probability distribution 
function (PDF) and the corresponding effective potential for several values of the coupling 
constant at the asymptotic future timelike infinity of de Sitter. For small field values the 
, effective potential has an approximately quadratic form, corresponding to a positive mass 

^^ ' term, such that the corresponding PDF is approximately Gaussian. However, the curvature 

o ■ 

■r^lj- ' of the effective potential shows qualitatively different (typically much softer) behavior on the 

o ■ 

• , coupling constant than that implied by the Starobinsky-Yokoyama procedure. For large field 

f^ ■ values, the effective potential as expected reduces to the tree level potential plus a positive 

cn ■ 

correction that only weakly (logarithmically) depends on the background field. Finally, we 

calculate the backreaction of fiuctuations on the background geometry and show that it is 

• 1—1 

^ \ positive. 

.^. 

I. INTRODUCTION 

The physics of de Sitter space is hard because there is no simple perturbative expansion parameter. In 
particular, massless scalars and gravitons exhibit so much particle production that their interactions cannot 
be perturbatively controlled. This breakdown of perturbative expansion is a serious obstacle to progress 
of our understanding of the physics of de Sitter space. The usual resummation techniques do not typically 
help. One such resummation is the self-consistent Hartree (mean field) approximation [1], which includes 
the one-loop resummation of daisy (and superdaisy) diagrams. When applied to an 0(N) symmetric real 
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scalar field on de Sitter space, this approximation scheme erroneously predicts that the 0(N) symmetry 
gets completely broken in the vacuum, and the (would-be) Goldstone bosons acquire a mass [1]. To get 
this result it is sufficient to assume the self-consistent Hartree and the de Sitter symmetry. A notable 
observation is that the mass of the Goldstones is strictly smaller than the mass of the condensate field. 
Recently, a more sophisticated resummation scheme based on a large N expansion in an 0{N) symmetric 
model has been utilized in Refs. [^, y| to study deep infrared correlators on de Sitter space. Furthermore, 
based on the Euclidean approach to de Sitter, progress has been also made towards understanding mass 
generation in an interacting scalar field theory with a quartic self-interaction and one real scalar field (the 
0(1) model) [J, l5|. Earlier resummation attempts include [6l4lll|. 

In this work we make use of a sophisticated resummation technique on de Sitter space known as 
stochastic inflation |l2l . Il3l |. Stochastic inflation provides a clever reorganization and resummation of 
perturbation theory in such a way that one obtains the correct leading order answers for infrared (super- 



12 



is that, while an interacting 



Hubble) field correlators. The main observation due to Starobinsky 
quantum field theory is essentially quantum on sub-Hubble length scales, it is classical on super-Hubble 
length scales if it couples non-conformally to gravity. Indeed, upon splitting the theory into short and 
long wavelenth modes, one finds that the dynamics of long wavelength modes is particularly simple: the 
modes with a super-Hubble wavelength exhibit overdamped dynamics, due to Universe's expansion gradient 
terms can be dropped and the only coupling between different modes comes from interactions. The coupling 
between the short and long wavelength modes can be modeled as a Markowian random force in the equation 
for the infrared modes. The resulting classical stochastic theory is particularly simple, and it can be shown 
to be equivalent to a Fokker-Planck equation for the single field probability distribution function (PDF) 
p = p{(f){x),t), which is of the form, 

dtP=^d4V'p) + ^dlp. (1) 

This equation is the stochastic equivalent of the von Neumann equation for the density operator in quantum 
field theory, and thus /)((/>, t) can be thought of as the classical limit of the density operator. 

Furthermore, it is known that a suitably adapted stochastic theory of inflation captures correctly fleld 
correlators in theories such as a self-interacting real scalar field theory [ij], Yukawa theory [l5| and scalar 



electrodynamics J16l4l8l]. Due to the complex interplay between the constraints and dynamical field com- 
ponents in gravity, no consistent stochastic approximation has been so far developed for theories that 
include dynamical gravity. This is a pity, since this leaves us with an incomplete understanding of the 
dynamics of quantum fields on de Sitter space. The task is further c omp licated by the fact that only a 
few perturbative results are known that include quantum gravity 19l428l|. The situation is much better 



as regards perturbative results for various field theories on de Sitter 29l438l|. Albeit useful, many of these 
results cannot be trusted at late times, when perturbative treatment breaks down. Therefore, one of the 
most burning unanswered questions of the physics of de Sitter space is: 

What is the effective field theory that governs dynamics and probability distribution of the infrared 
fields in interacting field theories on de Sitter space ? 

That question is particularly interesting in the context of eternal inflation, as there one expects that any 
time dependence can be viewed as a dependence on some physical scale. 

In order to make progress towards such an effective theory of inflation, note that - in the spirit of the 
renormalization group approach to effective field theories ~ the time dependence in ([1]) can be viewed as a 
(physical) scale dependence, where /i = /ioe~ (this dependence originates from the exponential expansion 
of physical scales in de Sitter space on super-Hubble length scales), such that dt — )• —Hfj^d^. With this 
observations Eq. ([1]) can be rewritten as a renormalization group (RG) equation, 

f^d.P + ^d^ {V'p) + ^dlp = . (2) 

This heuristic equation determines the PDF p = p{(f){x),p) as a function of the field (p and the scale /i, 
which should be interpreted as the physical scale on which the field varies. An important question is how 
the PDF changes if one integrates out fluctuations above the scale p. In the spirit of RG, one gets the 
effective PDF Pcs{4'i p) = Pfj.i'P) when fluctuations above that scale are integrated out. We shall postpone a 
concrete calculation of /J^((/>) oc exp( — |§^V"^((/>) ) for future work, and concentrate here on understanding 
the behavior of p^{(f>) and V^{(j)) in the limit when /i — )• 0. This PDF we refer to as the effective PDF, 
Peff, and it signifles the probability distribution of (deep infrared) flelds on the asymptotic timelike future 
inflnity of de Sitter space. The program we propose here is illustrated on the Carter-Penrose diagrams 
in figure [H The effective field theory (the effective potential and the corresponding effective probability 
distribution function) we discuss in this paper lives on the p = {t = c«) surface, and it is independent on 
the coordinates one uses. One can use this effective field theory to calculate the correlators of the fields that 
originate from (correlated) sub-Hubble vacuum fluctuations at some early time, and correspond to very 
large scale correlators of deep infrared fields at the asymptotic timelike future infinity of de Sitter. Albeit 
Eq. ([2]) is appealingly simple, it is not obtained by a rigorous derivation, and we shall address elsewhere 
the problem how to formally construct an effective field theory on de Sitter valid on some finite physical 
scale. 

The effective theory approach advocated here becomes even better motivated when one recalls the 
observed equivalence between the zero mode partition function Ze{4>oe) on Euclidean de Sitter space and 
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FIG. 1: The Carter-Penrose diagram of de Sitter space in flat coordinates used in this paper (left panel) and 
in global coordinates with positively curved spatial sections (right panel). The curves corresponding to constant 
physical scales ^^oo, ii — H,ii<H,fi<tiH and fi — are also shown. Even though flat coordinates cover only 
1/2 of de Sitter space (shaded), asymptotically the surface /^ — >■ (or equivalently, the future timelike infinity surface 
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on which t -> oo) is equal in both coordinates. We hence expect - and claim - that the effective theory that 



governs the field distribution on the surface fi = is independent on coordinates. 



the PDF in stochastic approach at asymptotically late times p{(j), t — ?■ oo) 
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where Ze denotes the Euclidean partition function, (JJOE is the zero field mode on Euclidean de Sitter 



space. 



60] Namely, while the Euclidean space zero mode is distributed according to Ze[4>oe\, the corre- 
sponding effective action, VEeflFi can be constructed by making use of a Legendre transformation. For a 
recent perturbative study of the Euclidean zero mode correlators see Ref. 5[. Inspired by this observation, 
an adaptation of this approach to stochastic inflation is the approach we advocate in this paper. 

In order to make our study specific, we shall consider an 0{N) symmetric scalar field theory, whose 
tree level action is of the form, 

N . N , / JV 
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where mo and Aq denote a bare field's mass and a bare quartic self-coupling, respectively, g^i, is the metric 
tensor, g^^"^ its inverse and o = det[g^u]- The metric signature we use is (—,+,+,..). Upon a standard 
renormalization procedure [l], tjiq — )■ m? = —p?, and Aq — >■ A > 0, where now m? and A are finite 
renormalized parameters. The resulting action can be used as a starting point for stochastic theory. When 
771^ < {p? > 0) the theory exhibits symmetry breaking in Minkowski space, which is what we assume 
throughout this work. In the case when A^ = 1, one obtains the action of a real scalar field, and we shall 



assume that the (renormahzed) action is of the form, 
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When studied on Minkowski space, the vacuum of the theory ([H) exhibits a scalar field condensate that 
breaks the 0{N) symmetry down to 0{N — 1) (in the case when A^ = 1, 0(0) means that the symmetry 
0(1) = Z2 is completely broken by the condensate). This spontaneous symmetry breaking (SSB) results 
in a non-trivial vacuum Ai = 0{N)/0{N — 1). Excitations along these vacuum directions are massless, 
and they are known as Goldstone bosons. For example, when A^ = 3, the vacuum corresponds to the 
two dimensional sphere, Ai = 0(3)/0(2) ~ S"^, such that there are two massless Goldstone bosons, 
corresponding to excitations along the two orthogonal directions on S^ . The vacuum of the theory is said 
to be trivial if (j){x) maps all of the physical space R^ into a point on field space, \\(f)\\ = 0o = ^y^iV/A 
(here A'^ = 3). If however (f){x) maps the 2-sphere ||x|| — )• 00 of the physical space onto a 2-sphere of 
the internal space, the vacuum is said to be topologically nontrivial (the second homotopy group of the 
vacuum manifold is non-trivial, 7r2(A^) = Z 39|). This configuration is known as the global monopole, 
and once formed it is (topologically) stable. Global monopoles are an example of topological defects, which 
have been studied in the 1980s and 1990s as an alternative to infiation that can seed Universe's structure 
formation. When confronted with modern cosmic microwave background observations j40||, these theories 



had to be abandoned. However, global defects have been invoked to drive infiation [41] or provide a possible 
explanation for dark energy [42| . 

The main purpose of this paper is to prove that no topological defects can survive in eternal infiation. 
In other words, deep infrared fields canno t see sy mmetry breaking and hence cannot form a condensate. 
For earlier work on this problem see Refs. 



43l445l|. The infrared effects from abundant particle production 
(in arbitrary number of space-time dimensions) in eternal infiation are so strong that they will eventually 
restore symmetry of any tree levelpotential of the form ([IHS]) and destroy any defects that might have 
formed on (sub-)Hubble scales J46l448l| . As an important side remark, we note that our results suggest that 
very deep infrared modes will exhibit an enhanced non-Gaussianity of a non-perturbative character. 

The paper is organized as follows. In section |II] we prove by explicit calculation that no condensate can 
form on de Sitter space in a theory (HD which at tree level exhibits a symmetry breaking {m? < 0). For 
pedagogical reasons, we present a separate proof for 0(1), 0(2) and the general 0{N) case. In section Hill 
we make use of the Legendre transform and show that the effective field theory at zero physical scale is 
convex for an arbitrary value of the background field. Again, we firstly study the 0(1) case, and then 
we prove convexity in the general 0{N) case. We also derive analytical approximations of the effective 
potential for small and large background field values. Finally, in section ITVl we summarize our main results, 



and point at future research directions. 

II. SYMMETRY RESTORATION 

In order to prove that the symmetry in an 0{N) model is restored, it is sufficient to show that the 
curvature of the effective potential governing deep infrared field fluctuations is positive at the origin. We 
shall prove this by reductio ad absurdwn. Namely, we shall show that the assumption that a condensate 
forms leads to a contradiction. We first consider a real scalar field, and then a complex (two component) 
scalar field and finally we discuss the general 0{N) symmetric case. 

A. A real scalar field 

Let us now consider a real scalar field, whose action is given in ^. We shall assume that the mass 
parameter /i^ = —m? > 0, such that the theory exhibits a spontaneous symmetry breaking in Minkowski 
space, and the field develops a condensate ip = ibi?i>o = i/xy^G/A in its trivial vacuum state. In addition there 
are (topologically stable) domain wall solutions in which the field asymptotically condenses to (/> = ib</>o for 
z — >• ibc« (where z is some spatial direction). Domain walls could have formed in the early Universe by the 



Kibble mechanism 39(], if at high temperatures the symmetry was restored by thermal field fiuctuations. 
In this case, as the Universe expands it cools down and the model undergoes a phase transition. 

We shall now consider the model ([5]) in infiationary (de Sitter) background. We shall argue that, 
while short scale field fluctuations can see a symmetry breaking potential, very long scale fiuctuations 
see necessarily a potential in which the symmetry is restored. Let us begin our analysis by varying the 
action ([5]) with respect to (p. On super-Hubble the derivative terms can be dropped and the equation of 
motion for the mean (background) field 0b simplifies to, 



y + -((^02)fi„ + -(/,2 

2 



0b = 0, (6) 



where the second term includes both the contribution from the fluctuations 6(j) = — 0b, that can be 
estimated by solving the Fokker-Planck equation ([1]) for stochastic inflation (while the sub-Hubble fiuctu- 
ations contribute mainly to renormalize the coupling parameters, the super-Hubble field fiuctuations are 
captured in ((50^)fin and can be estimated by stochastic infiation, (50^)fin = ((^0^)stoch)- 
Ignoring the 0b = solution (which is a local maximum), Eq. ([6]) implies, 

0^ = 02-3(502),toch>O. (7) 
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FIG. 2: The PDF (J)q x pi(0) as a function of <t>/(t>o for the Hubble scale modes. The three curves represent ^ = 0.1 
(solid blue double peaked curve), C = 1 (dashed red middle curve) and C, — IQ (dot-dashed green flattest curve), 
where C, = 47r^/i^/(Ai?'*) is the dimensionless parameter characterizing the inverse coupling strength. 



Assuming a stationary probability distribution p{(j),t) = p{(j)), Eq. ([T]) simplifies, and it is easily solved by 
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where C = (4vr^/i^)/(A/f'^) denotes a dimensionless 'inverse coupling' parameter. The second way of 
writing Vi in ([8]) is suggestive, as it indicates that - when written as a function of the rescaled field 
(p/cpo ~ all properties of the stationary p in (J17p can be parametrized in terms of just one parameter: the 
inverse coupling parameter C. Thus in the limit when ^ — )• oo the theory is weakly coupled, while in the 
limit when ^ — )• the theory becomes strongly coupled. While solution ([8]) has been known for quite a 
while [l3|, we give it a slightly different physical interpretation. For us Eq. ([8|) represents the PDF for 
the (stationary) modes that vary over approximately the Hubble scale. Albeit the PDFs that exhibit a 
nontrivial time dependence are interesting of their own right, for simplicity we do not study them here. 
Namely, Starobinsky and Yokoyama [13l | have shown that, after a sufficient amount of time, every initial 
state necessarily reduces to the PDF in ^, making the stationary distribution ([8|) an attractor. 

In figure [2] we show the initial PDF for these stationary (mainly) Hubble scale modes. We shall refer it 
to as the PDF for Hubble scale modes, and the corresponding V = V{(f), /x ~ H) is the effective (Hubble 
scale) potential. Of course, this potential still exhibits symmetry breaking, which can be seen from the 
characteristic double peak structure of the PDF in figure O which is more pronounced for strong couplings 



(C <C 1). Nevertheless, as we show below, the inclusion of these fluctuations into (^^ is sufficient to prevent 
background field condensation. Indeed, inserting 



into Eq. ([7]) results in 



)stoch = {(f> )stoch + 4>h 



i2 _ ^0 /i Q/ J,2 / i,2\ 



.b - ^ (1 - 3(</.V<Ao)stoch) > (9) 

_ 1 A , /3/4(C/2)+/-3/4(C/2) \ 
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where 



«-) = gf(*T7TTw(2) ("> 

is the series representation around z = for the modified Bessel's function of the first kind. When (jlOp is 
inserted into Eq. ([9]) one gets, 

which is impossible to satisfy, thus completing the proof. The proof is an immediate consequence of the 
observation that, if z^ > —1, Iu{z) is strictly positive for all z > 0. This is true because for i/ > —1 
all coefficients in the series (|lip for luiz) are positive when v > —1, and this condition is satisfied for 
all Bessel's functions occurring in (I12p . We have thus proved that, when the field fiuctuations evaluated 
within stochastic theory of inflation are taken account of, no phase transition in the 0(1) model ([5]) can 
occur, i.e. there can be no field condensate. 

B. The 0{N) model 

The proof in the 0{N) symmetric model dH is analogous to the proof in the 0(1) case. The main 
complication is the enlarged symmetry, and we shall devote some attention to explain the resulting differ- 
ences. In the 0{N) case the field (j) = {(pa) (o = 1) ••) N) is a A^-component vector, and when a mean field 
condensate is present, a suitable orthogonal rotation can bring the condensate direction into the direction 
of (/>i, such that 0b = (/'ib and (</>o)^ = (0i)o = Nfi^/X. The generalization of ([6]) to the 0{N) case is (c/. 
Ref. Ill) 

0^ = (0i)2 = 02 _ 3{{6<PifUoch - {N-l){{6<P2f)stoch > 0, (13) 



where 

((50i)2) = {(^,-^^f) = (cp^f + {(cp,f) = (^^^f + l-{($)^) (14) 

{{ScP2?) = {{^2?) = ^m')- (15) 

Because the scalar potential in (j3|) depends only on X]a=i('^a)^ ~ (.'P)'^i the linear term in ()14p does not 
contribute, and (((500^) = ((5(/'j)^) (Vi,i G {1, ..,N}). Inserting Eqs. (fTiHTSD into ([l3D results in 

0.^ = 1 (l- ^((^7^0) Vch)>0. (16) 

We shall now show that this inequality has no non-trivial solution, i.e. that no condensate can form for 
stationary solutions in stochastic inflation. 

But, before we consider the general case, it is instructive to consider the 0(2) model (complex scalar 
field), since in this case the results are particularly simple. We have 

(^ A^ 1 f ^""^A 1 12 VCe-i 4^2^4 

V2 = -^{^l + 4) + ^{^l + ^l? ^ ^^V2 = cWM'[WM'-2], (17) 

where erf (z) = {2/^/^T) J^ dte~^ denotes the error function and we normalized p by demanding 

fZo <^4>id4'2p = 1- 

In order to make progress on inequality (fT6l) . we need to calculate {4)'^)stoch- Making use of (fT7|l for 
N = 2 we immediately get, 

i4 /-oo 



j^4 /"OO 



e 



-C 



'°''^ ^[i;e,.f(VC)l j- ''^> 



where ip'^ = ('/'i+'/'D/'/'o ^^^ u = ip^ — 1. Upon inserting this result into ()16p with N = 2 yields 

2e-« 



,^ = -^ 



1 + 



>0. (19) 



/<[l + erf(VC)]_ 

Now, since erf(z) > (Vz > 0) Eq. (J19p cannot be satisfied for any C > 0, proving the impossibility of 
condensate formation for the 0(2) case. 

We are now ready to consider symmetry breaking in the general 0{N) symmetric case. In this case the 
PDF for a stationary state is the following generalization of Eqs. ([1] 
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(20) 



where r2(S'^ ^) = 2-k^i'^ /T{N/2) denotes the volume (surface area) of the unit A^ — 1 dimensional sphere. 
The two point function is then: 

= 0(5^-1)^/ du{u + l)^/h-<^\< 

VC r(f)xiFi(f;i;C)+2VCr(^)xiF,(^;|;C) ' ^ ^ 

where ip'^ = cp'^/cpQ and u = ip'^ — 1. Upon inserting (j2ip into (|16p we get 



a2 



/,2 



b = fF(C,iV) (22) 



where 



Ar + 2 2r(f + l)x,F,(f + l;i;C)+;^r(^)x,F,(^;l;C) 
^^' ^ A^ r(^)xiFi(^;i;C)+2VCr(^)xiFi(i^;|;C) ' ^ '^ 

In the appendix we show that -F(Ci N) < for any integer A^ > and C > 0, which completes the proof 
that there can be no mean field condensate in an 0{N) symmetric scalar field theory in de Sitter space. 



III. EFFECTIVE POTENTIAL 

While in section [TI] we show that inflationary fluctuations treated within the stochastic formalism are 
strong enough to prevent formation of a scalar condensate in any 0{N) symmetric scalar theory in de Sitter 
inflation, here we prove a much more powerful theorem. Namely, we show that the effective potential 
governing the PDF of very long wave length fluctuations in inflation is strictly convex (for arbitrary 
value of the field). This represents a general proof that deep infrared modes in de Sitter space see a 
symmetry restoring potential, and hence their PDF must be peaked at zero field value. Inspired by the 
observation that the partition function of the zero Euclidean mode equals the PDF for time independent 
field configurations ([3|) we shall use the effective action formalism to construct Vcs and the corresponding 
PeS- We first consider the 0(1) case, and then subsequently the general 0{N) symmetric model. 
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A. The 0(1) model 

Since p{(j)) yields a PDF for cj) that varies on the Hubble scale, it is reasonable to posit that pcS obtained 
by the conventional field theoretic technique of Legendre transform [49] will give a PDF for deep infrared 
(/U — )• 0) scalar field fluctuations 0b i whereby all higher energy fluctuations have been integrated out. Since 
any field that originates from sub-Hubble scale fluctuations rapidly redshifts during inflation, the effective 
potential that we derive here yields the field distribution at the asymptotically late times (t — )• oo), at the 
asymptotic future timelike infinity i~^ of de Sitter space, as illustrated in figure [TJ 

Adding a source current J(x) to V{(j)) one can define a partition function (at a point x) Z{J{x)) as 

Z{J) = e-'^^-^) = / #pi ((^)e-^'^ = (e^*) . (24) 



The effective potential V^g is then given as a Legendre transform of W{J) = — ln[Z(J)], 

|^T/eff(</)b) = H^(J) + J0b; ^^= ^\^/^ (25) 

and the corresponding PDF is then 

/5eff(</'b) = -^exp f-^^V;,ff((/>b)j ; Zes=d(l)bexpi-^j^Vesj. (26) 

In order to study convexity of I4flF((/>b) we shall make use of the equation of motion for (py, and its 
derivative, 

5<Ab ' 902 50b ■ ^^ 

V^ff is convex if d'^Ves/dcj)'^ > for all 0b- The inverse of this must also be positive, 

'9Veff(0b)\"' dHn[Z{J)] Z"{J) (Z'{J)^^ 



a2 



Now, from the definition of Z{J) in Eq. (j24p is follows 



> . (28) 



Z{J) = Y: ^ = {-'") - (29) 



n! 

n=0 



where 



i2n 



r(f + i) X iFi(f + 1; 1; C) + 2VCr(f + |) x iFi(§ + f ; |; C) 



^o" / .toch 7reC/2ct+i (/i/4(C/2) + /-i/4(C/2)) 



(30) 



and (02"'''^)stoch = (n = 0, 1,2, ..). Furthermore, we have 

Z'{J) = (0e-^'^) = Z(J)(0) J ; Z"( J) = (02e^<^) = ZiJ){cf>^)j , (31) 
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FIG. 3: Left panel: V^i as a function of 0b for C = 0.1,1,10 (solid blue innermost curve, solid red middle curve 
and solid green outermost curve). We also show the weak field analytic effective potential (15^1) (dotted curves), the 
analytic expression (|43l) for large (/)b (dashed curves), and finally the Starobinky-Yokoyama effective potential (|38)) 
(sparse green dashed innermost curve) just for ^ = 10 (in order not to overcrowd the plot). Right panel: Poff £^s 
function of 0b for C, = 0.1, 1, 10 (solid curves, the same color code as on the left panel) and the best Gaussian fits 
close to the origin (dotted curves). 

where (•)j denotes averaging with respect to pj{4>) = /o((/>)e'^'^, which is also positive. With this in mind 
we can rewrite the convexity condition ()28p as 

a2T- ' ' ^^ -1 



d'VM 



')J 



fj>0- 



(32) 



Now, since pj{(p) is positive definite, we can make use of the Cauchy-Schwarz theorem for probabihty 
theory, 



61| which for / = 1 and g = (p (see footnote [61]) reads 

2\ 



fj< 



(33) 



where the strong inequality follows because 1 and (p are linearly independent functions. This immediately 
implies that the (strong) convexity condition (j32p for Ves is satisfied, completing the proof. 

B. Analytical approximations for VeS 



In figure [3] we show numerical plots for Vcs and peS as a function of the deep infrared background field 
(/^bi for different values of the dimensionless inverse coupling parameter: C = 0.1, 1, 10. The figures do not 
represent exact V^s , but instead they are based on taking a finite number of terms in the sum (j29p involved 
in the definition of Z[J). However, we have checked that increasing the number of terms in the sum does 
not change the form of V^s to such an extend to be visible on the plots. It is apparent that, for small 
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values of 0b! Kff exhibits a simple quadratic dependence on 0b, implying an approximately Gaussian p^s 
(that can be characterized by a positive mass term), as can be seen on the right panel of figure [3l In 
the large field limit when 0^ ^> 0o the effective potential is, as expected, quartic, [{Sir'^ / {3H^)]Vcs ~ C'/'b- 
However there is no good smooth matching between the quadratic and quartic behavior, and around 
4>h ~ 4>0 the quadratic behavior turns quickly into a quartic behavior as one would expect from a Maxwell 
construction. Indeed, forcing a quadratic plus quartic fit onto VeS when C ^ 1 results in a fit with a double 
well structure. These observations then imply that weakly coupled theories exhibit large deviations from 
the behavior expected based on the Starobinsky-Yokoyama procedure, in the sense that small amplitude 
fluctuations of deep infrared fields are enhanced when compared with the Starobinsky-Yokoyama result 
discussed in section [TI] and below. 

In fact, it is worth making a detailed comparison of the mass implied by the effective action in figured 
and that implied by the Starobinsky-Yokoyama procedure [5|, ll3|] which, when generalized to a finite 
background field cpy^, yields. 



mSviM = -/^^ + 77((0 - 0b)^)stoch = /^^ 



k 


/*l^ 


\ ,1 


[ V 


\'Pi, 


stoch 'PO 



1 



(34) 



where (((/'/0o)^)stoch is given in Eq. pil|) . On the other hand, the (numerical) mass that measures the 
curvature of the effective potential in figure [3] is given by 



d^ 



M^ 



3W 



VesiM 



5(0b/0o)2 • 
Now, one can easily show that expressions pUj) and ([55]) are related as, 

m^(0b) 



^M'iM 



which - if the Starobinsky-Yokoyama prescription is correct - implies the following mass term. 



M 



SYl 



12C 



+ 



' stoch 



(35) 



(36) 



(37) 



The corresponding effective potential is then 



3H^ 



K 



effSYl 



.) = 6C 



1-'^ 






(38) 



/ stoch 

where Vq is an integration constant (a linear term in 0b is forbidden by symmetry) 

On the other hand, one can calculate the effective potential around 0b = by making use of Eqs. ([2 
25]) . Upon expanding e '^ in ()24p in powers of J0 and keeping the terms up to (J0)^) we get 



3^4 '^eflf - 


= 1 ^" +1 
2 (02) stoch 8 


^ 


l(0')stoch^ 


^t 


('^^)stoch 



+ 



(39) 



which yields for the mass around (j)]^ = 0, 
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^'(0) = 772V— ' (40) 

((P )stoch 

implying the following physical mass (j36p at ^^ = 

"^'(o) = rw^ — ■ (4^) 

There is a very intriguing connection of the result (j4ip . with the mean field (one loop) result for the 
mass [1[. Recall that in Z? = 4, for a field of the mass parameter m^, the mean field treatment gives 
{X/2)iA{x;x), where the coincident propagator is, iA{x;x) = (3/7^)/(87r^?n,^). Now, assume that the mass 
m is entirely created by stochastic fluctuations, thus m^ — ;■ (A/2)((^^)stoch; the one-loop mean field formula 
(with a stochastic flavor) then gives 7ti,^(0) = (3-fr^)/(87r^((/)^)stoch), which is precisely the result (HH). Of 
course, this is a very heuristic 'derivation', and it is very difficult to imagine that one could arrive at that 
result by any other but the rigorous method advocated in this work. 

At a first sight the mass (|40H4ip may appear bizarre, as it exhibits a very different dependence on the 
coupling parameter C oc 1/A from that implied by the Starobinsky-Yokoyama formula ()37p . At a second 
sight however (|^T]) is not surprising at all; it simply tells us that for small field amplitudes, 

Pes oc exp ( - — -^T^ ) ; ((/-b < M , 

which is precisely what one would expect from a stochastic theory in its Gaussian limit! (Albeit the theory 
is far from being Gaussian, for weak fields it is approximately Gaussian.) 

Let us now check whether our result ()40p agrees with the numerical findings in figure [3l In the 
weak coupling limit when ^ — )• oo, M^(0) — t- 1, while in the strong coupling limit (C — >• 0), M^(0) — t- 
[r(l/4)/r(3/4)]VC- For C = 0.1,1,10, Eq. (gOD predicts {0.72,1.20,1.03}, which agrees weh with the ap- 
proximate curves on the left panel of figure El These values are to be compared with the values implied by 
Eq. ([37]): M|y(0) ~ {1.3,6,77} for C = {0.1,1,10}. We have thus found that the Starobinsky-Yokoyama 
procedure typically predicts a much larger mass term at the origin than what is implied by the effective 
potential approach. 

In order to get a better feeling for the coupling dependence of the physical mass (j4T]l , consider the week 
and strong coupling regimes of Eq. ()4ip . 



2.o^A-^o Aff^ ^ A^ r(l/4) VAff^ ^ imJ^H^ (42) 

Thus, the strong coupling behavior is qualitatively the same as that of Starobinsky-Yokoyama [13] m?{0) ~ 
\/XH^ (albeit the dimensionless prefactors do not agree), but in the weak coupling regime the results are 
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qualitatively different. The latter result in Eq. (|42|) should not surprise us: in the strong coupling regime 
the curvature around the origin {—^J?) is negligible, and one recovers the expected result ?tt-^(0) ~ yXH^ 
(which is the loop counting parameter of the resummed perturbation theory). On the other hand, the weak 
coupling result in ()42p is quite surprising, in that it tells us that one almost recovers the naive dependence 
on the coupling constant, m?' ex \H^ (the additional factor H'^ / ^i^ is not easily explicable). What is 
interesting about that result is that, in the limit when A — t- and/or ^ — )■ cxd, small amplitude fluctuations 
at the timelike asymptotic infinity of de Sitter behave as massless field fluctuations, with the caveat that 
the amplitude of these fluctuations must satisfy \(t)\y\ < 0o = //y^G/A, which is the scale at which the (f>^ 
term in the effective potential kicks in. 

Next we consider the large field behavior of the effective potential ()25|) . In the Appendix we show that 
the asymptotic effective potential is of the form. 



^^^VesiM 


= c 


70by 2(^^Y 
\(j)oJ V0oV 




f ^AV,s{(t>b) 


|^AFefr(<Ab) 


C 1 /37rC0?\ / /C\ /C\\ 


or equivalently, 




v^siM = -^^l 


+ ^cPt + AV,s{M 


AV^siM « 1^ + 


167r2 
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[ H VH2J+^-H2JjJ 



^b > (po) 



(43) 



(0b»0o), (44) 

Firstly, the potential (|43H44p fits excellently the numerical Vcs for large (ph, as can be seen in the left panel 
of figure [3] (dotted curves). Secondly, the leading term {C{4'h/4'o)^) and the subleading term (— 2C((;^b/0o)^) 
are the same as in the tree level potential ([5]). However, there is also the contribution AV^fj that originates 
from integrating out fluctuations on all super-Hubble scales, and it consists of a logarithmic contribution 
((1/2) ln[((/)b/0o)^]) and a fleld independent contribution. These contributions are felt by large field exci- 
tations as a result of integrating out fiuctuations on all scales, and they introduce an average upward shift 
in the effective potential that grows weakly (logarithmically) with the field amplitude. This means that 
large amplitude, deep infrared fields in de Sitter are not distributed according to the Coleman- Weinberg 
effective potential 50]. A comparison with the Starobinsky-Yokoyama effective potential ()38p reveals that, 
in the large amplitude limit, it recovers the leading (quartic) contribution in the field correctly, but it fails 
at lower orders (quadratic, etc), and hence it does not represent a very good fit, as can be seen on the left 
panel in figure [3] (for (^ = 10 long dashed green). In order to get a better feeling for AT4ff in (|44|) notice 
that in the weak and strong coupling case it reduces to. 



Av^siM 



A^p Sfi^ , 3H^ ,„ fX(pl 



2A ^ 167r2 



In 



/"^ 



AT/ (^\ ^z^ ^^\ /^^VM_^ (Ar.^ 
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This means that the contribution from integrated fluctuations grows as A decreases. In fact there is an 
upper bound for AT4g which comes from the natural assumption, fJ,,(j)o ^ w-p; which then imphes that 
^^flf ^ "T-p/^, which is of the order of the Planck energy density. In the strong coupling regime however, 
AV^ff remains limited to AVcs ~ H'^. 

Next, we consider the backreaction on the geometry, for which we need the energy stored in the fluc- 
tuations. According to the Starobinsky-Yokoyama prescription, the energy density in stored in quantum 
fluctuations is 

(V')stoch = - Y('^^)stoch + Tj('A'*)stoch , (46) 

where ((/'^)stoch and ((^'^)stoch are given in Eqs. (lU\i and (fSO]) . For (" = {0.1,1,10} Eq. P6]) gives 
[87r^/(3-fr^)](y)stoch = {0.11,-0.58,-9.48}. On the other hand, we can numerically evaluate {Ves) and 
we get for C = {0.1,1,10}, [87rV(3i?^)](Kff)stoch = {0.39,0.34,0.23}, respectively. We thus see that the 
naive contribution from quantum stochastic fluctuations as given by the Starobinsky-Yokoyama procedure 
starts as a positive constant (= 1/4) when (^ — )• 0, but then it becomes negative as C increases, approaching 
[87r^/(3//^)](y)stoch — — 2C when C — >• oo. This is the energy density perceived by the fields that vary 
over the Hubble scale, but cannot be used for backreaction on the background geometry, for which we 
need the energy density perceived by very deep infrared modes, which is just (14ff)stoch- The contribution 
[Stt'^ / {3H'^)]{Vcs)stoch is (a) positive for all values of the coupling constants ( and (b) it decreases as C in- 
creases. This is so because, when A — )• and close to the origin where p is significantly different from zero, 
Vefi ^ and it is very fiat, and where Ves is large (0b ^ </>o), p shoots rapidly to zero. The energy (K!ff)stoch 
is to be added as a positive backreaction (87rGAr)(yoff)stoch to the background cosmological constant Aq, ef- 
fectively increasing the rate of Universe's expansion. Since [Stt'^ / {3H^)]{Vc{i) stoch is not greater than unity, 
the contribution to the cosmological constant is of the order AA = (87rGAr)(yofr)stoch = {^'^Gn)0{H'^), 
which implies AA/Aq = 0{H'^ /mi) and hence it is typically small. 

Next, we shall consider the entropy stored in fluctuations at a point x at asymptotic timelike future 
inflnity i^. Since stochastic formalism does not contain all of the information necessary to reconstruct the 
von Neumann entropy of the state 



Shannon (or Gibbs) entropy 
our problem yields. 



51 



54i |. we resort to the less fundamental, but simpler, concept of the 



511 ] ■ which is deflned as S'shannon = — (lii(p/io))stoch; which when adapted to 



5'shannon,SY = -^774 (^)stoch + In ( j , (47) 



which in the strong and weak coupling limit gives (for simplicity we set the scale pQ 



yoj 



5shannon,SY ^"^ In (4) + 4 + In (^Tf^^J ^ ^Shannon.SY '^'^ " 3 1" {' 
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The fact that in the weakly coupled regime 5'shannon,SY < should not be of a concern, since S'shannon.SY 
in (j47l) is defined up to a constant determined by a mass scale ^q, and ^q can be chosen such to keep 
'S'shannon.SY > 0. In particular, for ( = {0.1, 1, 10} one gets, 5shannon,SY = {1-52, 1.09, -0.033}. This is to 
be compared with the Shannon entropy implied by the effective theory, Sshannon,eff = —{^^{PcSl-i'o))stoch, 
which yields S'shannon,eflf = {1.375,0.983,0.812} for ( = {0.1,1,10} and fiQ = cpo, respectively. Of course, 
•Sshannon.cflf > for arbitrary coupling constant provided fio is chosen large enough {fiQ = (J)q suffices). This 
is so simply because Vcs > (V(/)b> C > 0)- Again, we have seen that the Shannon entropy of the effective 
field theory gives more reasonable results than the corresponding quantity calculated by the Starobinsky- 
Yokoyama procedure. 

Finally, we give one last cursory look at the V^s shown in figure [3l and observe that in the weak coupling 
limit (see the (" = 10 curves on the left panel of figure [3j) the effective potential Ves become fiatter around 
the origin {(pi, < (po) as C becomes larger, and it becomes steeper for large values of the field ((/)b > (po), such 
that there is a sudden change in the effective potential curvature at 0b ~ 4'o- These features are reminiscent 
of a Maxwell construction (see e.g. Ref. [49y), which is an approximate procedure for constructing the free 
energy/effective potential close to a critical point (where fluctuations become massless and long range 
correlations develop). Indeed, in the weak coupling regime one can approximate the effective potential 
by the following Maxwell-like construction: use Eq. ()39p for i?i>b < <Po a-nd Eqs. (I43I444I) for (p^, > (pQ, and 
this will give a reasonable approximation to the true effective potential, that exhibits a discontinuity in 
curvature at (/)b ^ (/>o. 

C. The 0{N) model 

In analogy to the convexity condition in the 0(1) model ()27p . in this more general case the effective 
potential Vcs{(j) ) {(p = {(Pi,i = 1, ..,-/V)) will be convex provided the corresponding Hessian matrix 



b'^^dcp) d(t>) 



(48) 



is strictly positive (i.e. all of its eigenvalues are strictly positive), where J = {Ji) is the source current 
{pf{(p) = p{<p)e '^). Taking the inverse of the matrix (|48p (which exists provided its determinant does not 
vanish) we get (c/. Eqs. (I2SH3D 

i^~^^) =(^^^^h-i<t'^).ri<t>^).j- (49) 

The simplest way of finding the eigenvalues of the inverse Hessian matrix (j49p is to make use of the 0{N) 
symmetry, which implies that there exists an orthogonal transformation R {R • R^ = I = R^ • R) which 
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brings the source current J to the first component J' = R ■ J , where J' = (J, 0, ..,0) and J = \\J \\. Now, 
making use of the invariance J • (/) = (J') • <j)\ <j)' = R • (p, and of the fact that the Jacobian of the 
transformation </)—)•(/)' equals unity, we conclude that the partition function Z{J) = Z{J') (this just means 
that Z{J) is a function of J^, which is invariant under orthogonal transformations). Consequently, we 
can act with R^ from the left and with R from the right on (j49p to obtain a rotated Hessian matrix, 
{'P'i^'j) J' ~ {'P'i) pi^'j) 7" which is in its diagonal form! In order to see that, we shall now write the positivity 
condition for the nonvanishing elements of the rotated inverse Hessian matrix, 

where in the last inequality we used the fact that {(k'i) v = for j = 2,3, .., A^, which is a consequence of 
the fact that J'^ = (i = 2, 3, .., N). Now the first inequality is proved identically as in the 0(1) case, while 
the second line inequality is trivially true. This completes the proof that V^g is convex, implying that deep 
infrared scalar fields <f)^ on de Sitter cannot undergo phase transitions. 

Needless to say, one can perform numerical and analytical analysis of the resulting T4g {4'^) by using 
analogous methods as in sections IHI Al and IHIBI 

IV. DISCUSSION 

We use the formalism of Starobinsky's stochastic inflation to show in section |II] that there can be 
no scalar field condensate formation in an 0{N) symmetric scalar field theory. For simplicity we have 
separately considered the 0(1) model (a real scalar field) and the 0(2) model (a complex scalar field). 
Nevertheless, as figure [2] clearly indicates, the probability distribution function for stochastic field that 
varies over the Hubble scale shows a double-peak structure, implying that first order transitions can exist 
during infiation, such that one can foresee formation of topological defects on sub-Hubble scales that persist 
on super-Hubble scales, at least for a while. 

In section IIIII we borrow the concept of effective field theories to study the probability distribution of 
the deep infrared fields in de Sitter space, that 'live' on the asymptotic future timelike infinity i'^ of de 
Sitter. We prove that the effective potential that governs the probability distribution of deep infrared fields 
of an 0{N) scalar field theory must be strictly convex. Our proof of convexity holds for arbitrary integer 
A^ and for arbitrary values of the background field. Our proof applies for an arbitrary quartic (Hubble 
scale renormalized) potential, and we strongly suspect that, by making use of analogous methods, one can 
construct a proof of convexity for an arbitrary 0{N) symmetric local scalar field theory. An important 



19 

consequence of this result is the inipossibihty of field condensation at the asymptotic future timelike infinity 
i"*" of de Sitter space illustrated in figured! The physical reason for this symmetry restoration is the strong 
infrared vacuum fluctuations on de Sitter space. 

It is worth mentioning that the effective potential can be fully characterized in terms of just one physical 
dimensionless parameter (^ defined in (I20p . which characterizes the inverse coupling strength. Furthermore, 
while for small background field values the effective potential is Gaussian (corresponding to a positive mass 
parameter [5|, ll3|]) for large field values it reduces, as expected, to the tree level quartic potential plus a 
correction that logarithmically depends on the background field. On the other hand, as our analysis in 
section IlII Bl shows, the curvature of V^s is much softer than the one implied by the Starobinsky-Yokoyama 
approach. 

In section IIIIBI we also discuss the backreaction on the background geometry, and show that it is always 
positive and suppressed independently on the choice of parameters. The backreaction is typically of the 
order AA/A = 0(A/mp), which is to be contrasted with the result obtained by the Starobinsky-Yokoyama 
procedure, which can be either positive or negative (in the limit when (^ — t- do, (y)stoch — ^ ~2^). 

Since astronomical observations of cosmic microwave background anisotropics and Universe's large scale 
structure measure some fixed finite physical scales, it would be of interest to extend the analysis in this 
paper to study a real scalar field probability distribution, and the corresponding field correlators, on some 
finite physical scale fi <^ H. In order to properly study these, one would have to develop a full effective 
theory for the probability distribution pcs{(j), fj-,t) on some finite scale //, for which a first - but possibly 



naive 
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a, ISSj] - guess is given by Eq. ([2]). Since Pes{4'-, /^i t) is the classical equivalent of the density operator 
of quantum field theory, peff{<j), fi,t) must contain complete information about the theory, that includes 
(equal time) field correlators, which are of relevance for cosmology. 

Finally, it would be of interest to extend our analysis to other stochastic theories whose evolution is of 
the Langevin type. Interesting examples include the infrared dynamics of thermal field theories such as 



non-Abelian gauge theories 
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58|. 



Appendix 
0(N) case 

Here we show that the function -F(C; A^) in Eq. (I23p is strictly negative for any positive integer A^ and 
real C > 0. Let us begin by recalling the series around z = of the confiuent hyper geometric function 
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(which has infinite radius of convergence): 



^. , N Ar(a + n)r(6)z" 
^T{a)T{b + n)nl 



(51) 



where T{z) denotes the gamma function, T{z + 1) = zT{z). It is now convenient to rewrite NF{(,N) 
A{C,N)/B{C,N), where: 
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(52) 
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(53) 



Let us first consider i?(C, N) in Eq. (j52p . In order to show that B > for any C > and for any positive 
integer A^, observe first that the confluent hyper geometric function iFi{a;b; z) >Oifa>0, 6>0 and 
z > 0. This is so simply because all of the coefficients in the series (|5ip contain gamma functions of positive 
arguments, which are strictly positive. Now, because Eq. (j52p contains confluent hypergeometric functions 
with positive indices a,b and of a positive argument ^ > 0, it follows that B{(^,N) > for any C > ^^^ 
any positive integer A^. 

To complete the proof we need to show that the function A in (I53p is strictly negative. In order to show 
that, observe that one can combine the first with the third and the second with the fourth term, to obtain. 
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(54) 



which is strictly negative since all three terms are negative for arbitrary C > ^-nd positive integer A^, 
completing the proof. Notice that the first term in (I54p originates from the first term in the series of the 
last hypergeometric function in (|53p . and that all gamma functions appearing in (|54p are strictly positive. 



Asymptotic effective potential 



In this appendix we present a derivation of the effective potential (j43p for large field values (ph ^ 4'0- 
In doing so we shall make use of the formulae (|24H26p and (l29H30p in the main text. Since the mapping 
J I— 7- 0b is a monotonically increasing function, when 0b is large so is J. When summing the series (I29p . 
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observe that in the hmit when J is large, terms with a large n dominate the sum, and hence we need a large 
index expansion of the confluent hypergeometric function, which can be, for example, found in Eqs. (9.228) 
and (9.220.3) of Ref. [59|, 



laKoo r(/3)e-/2z^ / ///3 

ihi[a]p]z) ~ ■;= X cos I 2a/ I — — a |z 



from which we infer, 

Inserting this into ([29H30|) yields 
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r(2n + l) 

where we have used a rescaled (dimensionless) current j = (J)qJ (similarly, we shall be using a dimensionless 
background field (^b = (ph/fpo)- The sum (I57p is hard to evaluate. However, by noticing that in the limit 
of large j, the sum is dominated by large n's, we can replace the sum by an integral, Yl'^=o ~^ /o ^"'' ^^^ 
make a saddle point approximation to the integral. This means that we can use a large n expansion of 
'^('^)C;j)) and approximate it by its expansion around the stationary point no, where dco/dn = 0, 



w(ra,C,j) ~ wo(C,j) + oWo(C,i)(^-^o)^ 



(58) 



where (in the limit when n — )• oo), 
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(59) 



The resulting Gaussian integral is convergent if u}q{C, J) > 0, which is indeed the case, and it evaluates to. 



Z{j) = e 



-^(J) Rs e"'^" 




1- -erfc(^/u;o72 x uq 



(60) 



where the complement error function erfc(a:;) = 1 — erf (x) can be neglected when the argument is large 
because, in the large argument limit, erfc(x) ~ e~^' /[-^/ttx]. The standard expression now yields the 
following approximate expression for 93b; 
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(61) 
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which when inverted gives, 

J « AC^l [l - i + 7X74 1 • (62) 



Inserting this into Eq. ()25p and making use of (j60p yields the asymptotic effective potential 



^TT' 



Kfr(0b)«C(v^^b-2c^^ + ^)+^ln(^)+ln(/.(0+/_.(0), (.^^ » 1) , (63) 



3i?4 
which is used in section IIIIBl see Eq. (jl3|) . 
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